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ABSTRACT: 


A restoration of the 37 „ 21, 18, 10.7, and 6.6 GHz satellite imagery from 
the Scanning Multichannel Microwave Radiometer (SMMR) aboard Nimbus-7 to 
22.2 km resolution is attempted using a deconvolution method based upon 
nonlinear programming. The images are deconvolved with and without the aid of 
prescribed constraints, which force the processed image to abide by partial 
a priori knowledge of the high-resolution result. The restored microwave 
imagery may be utilized to examine the distribution of precipitating liquid 
wter in maritime rain systems. 

1. INTRODUCTION 

Knowledge of maritime precipitation amounts is useful in climatological 
energy budget studies, as well as in the diagnosis of latent heating processes 
in individual storms, such as tropical cylcones. In addition, remotely-sensed 
rainfall distributions that are retrieved in real time can be incorporated 
into forecast model initialization schemes. 

The sensitivity of upwelling microwave brightness temperatures at the 
^canning Multichannel Microwave Radiometer frequencies (37, 21, 18, 10.7, and 
6/6 GHz) to variations in precipitation liquid water content of maritime 
clouds has been demonstrated both theoretically (Olson, et al., 1982) and by 
statistical correlation of measured radiances with available radar data 
(Spencer, et al., 1983). This sensitivity provides a means to remotely 
determine precipitation liquid water content in storms over the ocean, where 
measurements from radar or rain guage networks are unavailable. 

Due to the fixed aperture size of the SMMR, the spatial resolution of the 
instrument at the earth's surface diminishes in inverse proportion to the 
channel frequency. Thus at 37 GHz the footprint dimension is 18 km x 27 km, 


f 


compared to a maximum footprint dimension of 95 km x 148 km at 6.6 GHz. It 
has been demonstrated theoretically by Olson, et al. (1982) that the lower 
frequency SMMR channels are more sensitive to changes in cloud precipitation 
liquid water content than the 37 GHz channels at moderate to high rainfall 

rates, if all measurements are made at ideal resolution. Unfortunately, much 

of the information contained In the lower frequency channels is degraded due 
to the inadequate spatial resolution of these channels. This kind of image 

degradation is especially serious in coastal regions, where blurring of a 

bright coastal land mass can obscure rain features over the adjacent ocean. 

In the present study, digital image processing techniques are applied to 
the SMMR imagery for the purpose of restoring all channels to a common 22.2 km 
resolution. 

2. SOLUTION OF THE DECONVOLUTION PROBLEM 

The effect of finite antenna aperture size on sensor resolution can be 
modeled as an antenna blurring function, , applied to an ideal resolution 
image f. The antenna blurring function is estimated in the next section for 
each of the SMMR channels. Given the effect of additive noise e, the 
resulting low-resolution, observed image g is expressed 

g = fi * f + e . (1) 

Here, the symbol * denotes the two-dimensional convolution operation. 

The relationship (1) actually represents a system of equations governing 
the unknown f in the discretely-sampled image domain. If the observed image 
g is interpolated to an I x J cartesian grid, and H is a square matrix with 
dimension N, then there exists for each image grid element g. , a governing 

' J 

equation 
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s = i - C(N-l)/2] + k-1 


t = j - C(N-l)/2] + £-1 


where the unknown f is an (I+N-l) * (J+N-1) cartesian image grid. 

Both the estimation of the antenna blurring function ft and the deconvolu- 
tlon of g to estimte f arc. accomplished using the nonlinear programming 
method described by Yeh, et al. (1984). Basically, the method considered by 
Yeh optimizes a given function (f or H in this study) by minimizing the mean- 
square error 
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subject to tne constraining relation 


P(f, ft) = 0 . (4) 

The function P(f, H) expresses the difference between some property of the 
current estimate of f (or ft) and that of the desired result, based upon 
physical laws or other a priori information. Specific constraint functions 
employed in the estimation of the antenna blurring function and the 
restoration of SMMR imagery are defined in sections 3 and 4b, respectively. 

A nonlinear programming method referred to as the penalty function method 
by Bazaraa and Shetty (1979) is utilized to solve the constrained optimization 
problem (Eqs. 3 and 4 above). In the penalty function approach, the 

constrained optimization problem is treated as a sequence of unconstrained 
problems. For each step of the sequence the combination 
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Is minimized using a standard technique to optimize f (or H), The conjugate 
gradient method of Fletcher and Reeves (1964), which allows faster conver- 
gence than the steepest descent method for differentiable functions such 
as i|»(f,H,A), was employed to optimiate f (or H) In this study. 

In the initial step, the parameter A is set to a small value, such that 
the norm of the penalty A • P(f,fi) is negligibly small with respect to the 
norm of 4i(f,fi,A). A value of A = 10' 5 was determined to be sufficiently small 
for applications in this study. Since the Initial weighting of the constraint 
function P(f,H) is relatively small, subsequent minimization of if>(f,H,A) using 
the conjugate gradient method leads to an optimal f (or fi) which is close to 
the unconstrained result. 

On each succeeding step the parameter A is increased by a factor of ten 
following by minimization of i|i(f,i$,A) by the conjugate gradient method. In 
this way, the influence of the constraint is gradually introduced. 

It is demonstrated in Bazaraa and Shetty (1979) that under certain 

general conditions the norm of the penalty A»P(f,H) tends to zero as A tends 

4 > + 

to infinity by this procedure, and the norm of <Mf,H,A) meanwhile approaches 
a local minimum. In this study the sequence of steps was terminated 

3 

with A = 10 on the last step, since it was observed, in general, that the 
norm of the penalty became negligibly small for this step, and the norm 

of i|»(f,H,A) converged to a minimum value. 

Since the image deconvolution problem is generally underdetermined (see 
Section 4a), it is possible that the parameter space of <f(f,H,A) has numerous 

local minima, and that the solution f obtained by local minimization 

•f 

of t|»(f ,H,A) may differ from the ideal-resolution image f. The imposition of 
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restrictive constraints will reduce the number of possible solutions. Also, t 
choosing an initial guess for f which is close to the expected solution (e.g. 
the original image g), spurious solutions may be avoided. 

If desired the unconstrained solution may be obtained simply by setting A 
equal to zero initially, 

3. ESTIMATION OF THE ANTENNA BLURRING FUNCTION 

In order to determine the antenna blurring function H, cloud-free SMMR 
swaths of the Florida coastline were first selected. Based upon physical 
models (Olson, et al., 1982) one would expect that the microwave radiances 
emitted from land and ocean areas to differ by at least 50 K. In addition, 

i 

the absence of clouds suggests that in each region the microwve brightness 
should be almost spatially uniform. Only minor brightness variations would be 
introduced by inhomogeneities in atmospheric water vapor and, in land regions, 
differences in soil moisture. At ideal resolution, therefore, the land/ocean 
interface represents nearly a step-function discontinuity in brightness 
temperature. The effect of antenna blurring is to smooth this discontinuity 
in the observed imagery. 

Using the cloud-free SMMR imagery, an approximate ideal-resolution image 
(f in Eq. 1) was created for each channel. The corresponding low-resolution 
brightness temperature data in scan format were first interpolated to a local 
6° latitude x 6° longitude cartesian grid with 22.2 km grid spacing (5 grid 
point per degree latitutde) using the Laplacian smoothing spline fitting 
procedure of Wendelberger (1981). Land and ocean brightness temperatures, 
assumed to be spatially uniform at high resolution, were then estimated from 
brightness temperature values sampled far from the blurred coastline. 
Coastline geography from USGS survey maps and the 37 GHz imagery were utilized 
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to locate the land/ocean brightness discontinuity at grid resolution. An 
ideal-resolution brightness temperature image f for the 37, 21 18, 10,7 and 
6.6 GH channels in both parallel and perpendicular polarizations was 
estimated in this way. 

4 4 

Given f and the observed image g in each channel, Eq. (1) was solved 
for using the method of Yeh, et al . (1984), subject to the condition 
that ft should be normalized* 


♦ N 2 
P (*) " (( I H ) - 1) Z 
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Each blurring function was defined on an N x N moveable grid, or 
convolution template, where the dimension N was adjusted to accommodate 
blurring operators of various sizes: for example, the 37 GHz channels, with a 
footprint size of 18 km x 27 km, required a 5 x 5 template (110 km x 110 km 
22,2 km grid spacing), whereas the 6.6 GHz channels required a 15 x 15 
template (330 km x 330 km). The convolution template for each channel was 
made large enough to ensure that the magnitude of the blurring funciton 
outside the template boundaries was negligibly small. The template dimension 
was therefore chosen such that the blurring function value at the template 
edge was less than 10% of the central maximum value. A schematic depiction of 
the 37n GHz convolution template is presented in Fig. 1. 

Separate blurring functions for day and nighttime SMMR overpasses were 
determined, since the scan line orientation relative to the image grid axes 
differs for ascending (daytime) and descending (nighttime) swaths. The 
initial guess blurring function for each channel was a Gaussian field with a 
dispersion width (2o) roughly equal to the published footprint dimension (ref. 
Gloersen and Hardis, 1978). 
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In Pig. 2, the estimated blurring functions for the daytime SMMR overpass 

t 

In the perpendicular polarization at c7, 18, and 6.6 GHz are presented. It 
should be noted that the antenna blurring functions derived here Include the 
average effects of varying scan geometry within the swath and sensor 
integration time. The resulting functions are therefore smoother than those 
that would be expected from aperture diffraction theory (see Jackson, 1975, 
p. 441 ff.). 

The figure dramatically illustrates the large difference in the magnitude 
of blurring between the highest and lowest frequency SMMR channels. Because 
the observed Imagery results from a convolution of the Ideal-resolution 
imagery with the blurring function for each channel, the sharply peaked 37 GHz 
antenna function yields observed imagery with spatial resolution only sligntly 
worse than the grid resolution of 22.2 km. Conversely, the 6.6 GHz blurring 
function is nearly flat, which leads to observed imagery with relatively poor 
spatial resolution. 

4. APPLICATION OF THE DECONVOLUTION METHOD TO SMMR IMAGERY 
4a. The Model Constraint 

Given the estimated blurring function H, the system of equations which 
must be solved for the high-resolution image f is necessarily underdetermined 
as suggested by Eq. (2). This arises from the fact that the convolution 

4 , 

template is of finite size. The convolution by H outside the boundaries of 
the observed image g. The number of equations in the system is equal only to 
the number of grid elements in ?, which always has a greater dimension, This 
situation is illustrated in Fig. 1. The larger the blurring template H, the 
larger the border region in f, and the more degrees of freedom the system 
acquires. 
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The degrees of freedom In the system can be reduced If certain partial 
knowledge of the solution can be determined a priori . An analytical model for 
microwave brightness temperature, upwelllng from a precipitation-free 
atmosphere over the ocean has been developed to provide such prior 
knowledge. The model is based upon a simplified version of the radiative 
transfer solution for the brightness upwelllng from a plane-parallel, 
isothermal, absorbing atmosphere. The brightness temperature at viewing 
angle 6(j = cos"^ at frequency v and polarization p is: 

TI V u *>V T sfc> ■ o-'V-Wv < 7 > 


+ "*atm v * ^ 1+r vp fi vp n v^ 


+ WvpVv 


Here, T sfc is the ocean skin temperature, n v is the total transmittance of 
the atmosphere and T bb is the cosmic background radiance {2.7 K). The sea 
surface reflectivity r vp is given by the Wentz (1983) parameterization. 


r vp(T S fc- u * - - V* • (8 > 

where e vp , the calm sea surface emissivity, and M , and emperical 
constant, are both given by Wentz, and U* is the surface friction velocity. 
The variable n vp , parameterized 



1 + u *U 

up 
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(9) 
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with <*> vp an empirical constant, Is another factor provided by Wentz to 
account for the change In sea surface reflectivity due to wind-roughening. 


Finally, the effective atmospheric emission temperature 
parameterized 

* K «T . 
atm^ v sf c 



Is 


( 10 ) 


The emperical constant K v were determined by fitting the model to clear- 
atmosphere SMMR brightness temperature data from the Gulf of Mexico. Values 
of K y , listed In Table 1, display a slight dependence on clear-atmosphere 
optical thickness, which is a maximum at the 22.235 GHz water vapor resonacne 
line. The atmospheric transmittance is further decomposed as 


"v ■ "dr v ,n cl<J v , 


( 11 ) 


where n c -| r is the transmittance is clear air, and 


n cl<j v * (n cld 37 > 


( 12 ) 


is an approximation for the transmittance through a layer of non-precipitating 
cloud droplets. Equation (12), with constants y v listed in Table 1, is 
based upon Rayleigh theory which accurately describes absorption by droplets 
with small size parameters x = 2trac/v. Here, a Is the droplet radius, c is 
the speed of light in a vacuum, and v is the frequency of radiation. A size 
parameter upper limit of .02 is obtained for a non-precipitating droplet 
radius of 100 urn at 37 GHz. The Rayleigh criterion is therefore well 

satisfied. 

Utilizing Eqs. (7-12), it was possible to estimate the 21, 18, 10.7, and 
6.6 GHz brightness temperatures in the rain-free area of each scene at 37 GHz 
resolution by the following method. The SMMR images in all channels were 
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first gridded at 22.2 km resolution as described In the preceding section. 
The rain-free portion of a given scene was then determined by thresholding a 
histogram of the gridded TB 37| - TB 37l values, which are Inversely 
proportional to the atmospheric optical thickness. The TB 37) - TB 37i 
distribution was found to be blmodal In all cases, with a distinct modal peak 
comprised of clear or nearly clear Image grid elements over the ocean in 
addition to a broader, less well-defined peak generated by grid elements In 
rain regions or those over land. A double Gaussian fit to this distribution 
was performed. Rain-free grid elements were defined as those for which the 37 
GHz polarization difference fell within one standard deviation of the "clear- 
air" modal peak and also those grid elements with greater TB 37(1 - TB 37jl . 

The observed parallel and perpendicular 37 GHz brightness temperatures 
and Eqs. (7-12) were next employed to determine n 37 and U* for each image 
grid element in the rain-free area. T g ^ , the only other unknown In these 
formulae, was derived from auxiliary ocean buoy data. In an operational 
framework, T $ f C might be approximated by the climatic mean temperature or 
estimated from satellite infra-red measurements. The inversion of n 37 and 
U* was performed numerically using Newton's method In two dimensions. 

Image grid elements in the rain-free area which contained the maximum 

values of n 37 were assumed to be cloud-free. Given U* for these elements and 

observed brightness temperatures in the other SMMR channels, estimates 

of n i were made using Eqs. (7-12). Although it may be argued that antenna 
v 

blurring alters the brightness temperatures in the lower frequency channels, 
the clear-air transmittances locally calibrated in this manner were found to 
be more accurate than those derived from empirical relationships. 

Having determined Jl 37 and n c ^ r for a given image, n y could be computed 
for any grid element in the rain-free area using the inverted value n 37 and 
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Eqs. (11) and (12) The 21, 18, 10.7, and 6.6 GHz brightness temperature at 
37 GHz resolution were then generated using n y , U*» and Eq, (7) in the rain- 
free region. 

The method so described Is essentially an extrapolation procedure, where 
variations in microwave brightness are mainly attributed to variations In non- 
precipitating cloud amount. In Section 4b, the 37 GHz resolution data Is 
utilized to constrain the image deconvolution system. 

It should be noted that although the histogram method, described above, 
discriminates reasonably well between grid elements in raining and rain-free 
regions, some ml selassif* cation of grid elements does occur. To Illustrate 
this point, histograms of grid elements In both raining and rain-free regions 
over the Gulf of Mexico and Florida coast at 0540Z, 7/17/79 are presented in 
Fig. 3. Rain areas were identified using grldded WSR-57 radar data from the 
station at Apalichlcola, Florida. Details of the radar gridding procedure are 
discussed In Section %. The single standare deviation threshold Is indicated 
by the vertical line at the 37 GHz polarization value of 54.1 K. 

The single standare deviation threshold was designed to produce a fairly 
conservative constraint. If this threshold is used then the constraint model 
is applied to grid elements which are almost certainly rain-free, despite some 
mlsclassification of rain-free grid elements. Nevertheless, grid elements in 
rain regions are sometimes classified as rain-free, as it is Illustrated in 
the figure. This latter type of misclassification is more serious because the 
constraint model was formulated under the assumption of a rain-free 
atmosphere, and might otherwise produce erroneous constraint brightness 
temperatures. In the future, correct classification of grid elements could be 
achieved using additional high resolution satellite data, such as that 
provided by visible and infra-red sensors. 
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4b. Restoration of SMMR Imagery 

Eleven 6° x 6° test areas over the Gulf of Mexico and eastern Florida 
coasts were chosen to evaluate the Image restoration method. The areas 
selected (a) contained significant maritime rainfall, (b) were located within 
range of operational WSR-57 coastal radar, and (c) were located within the 
swath boundaries of the Nimbus-7 SMMR during overpasses in 1979. The 

brightness temperature data, In scan format, were first interpolated to 
latitude/longitude cartesian grids at 22.2 km spacing as described In 

Section 3. 

Based upon the observed Imagery and estimated antenna blurring function, 
Eq. (1) was solved for the high-resolution Image f in each channel, using 
nonlinear programming method described by Yeh, et al , (1984). The model 
cons.ivlnt, presented In Section 4a, was Incorporated into the constraint 
functions P(f), 

p( *> ■ w 2 (13) 

* »J 

summation over rain-free area 

where TB^j is the model -generated brightness temperature. 

Two variations of the restoration method were considered. In Case A, the 
restored image solution was left unconstrained (A=0), and was initialized 
using the observed image throughout the domain of g. Thus no prior knowledge 
of the solution was utilized in Case A. Case B restorations were fully 
constrained in the rain-free portion of the image by allowing A in Eq. (5) 
to increase in steps as the iteration proceeded. The initial guess in this 
case was the observed low-resolution image except in the rain-free portion, 
where the model constraint was imposed. In both Case A and Case B the initial 
guess was extended to the border region of f by repeating the edge grid 
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elements In the domain of g outward. Note that since the model constraint 

was derived from the 37 GHz data to be applied only to the other SMMR 
channels, Case B was omitted for the 37 GHz Imagery. 

Examples of restored Imagery from the two deconvolution cases are 
presented in Fig. 4, where the original Imagery has been Included for 
comparison. The scene In the figure covers nearly the full extent of 
Hurricane David at 1644Z on 9/3/79, as It made its landfall on the eastern 
Florida coast at out 90 km south of Datona Beach. 

Despite the underdetermined nature of the governing system (Eq. 2), 
satisfactory results were obtained at 37, 21, 18, and even 10.7 GHz when the 
system was left unconstrained (Case A). The missing features and artifacts 
evident In the Case B restorations can mainly be attributed to error in the 
spatial constraint, arising from both misclassification of "rain-free" image 
elements and from error in the constraint brightness temperatures themselves. 

Other types of constraints on the restored-resolution satellite imagery 
were considered by the authors. These included imposed maximum and minimum 
bounds on brightness temperature, a constraint on the total image brightness 
to remain constant, and constraints on the intra- and inter-channel grid 
element brightness temperature cross products. The later constraints were 
specifically designed to preserve the mean auto- and cross-correlation 
statistics of the grid element brightness levels. Unfortunately, none of 
these constraints resulted in significant improvement of the restored image 
quality. It should be noted that minimization of ^(f,H,A) (Eq. 3) results 
in preservation of the expected value of the grid element brightness 
temperatures . 


5. COMPARISON OF RESTORED MICROWAVE IMAGERY WITH WSR-57 RADAR 


For each of the eleven gridded SMMR scenes a dlgnitlzed radar map in the 
same grid projection was produced. Radar contours from coastal t.'3R-57 radar 
were transferred to a grid overlay, where each 22.2 km grid box of the overlay 
was centered on a corresponding SMMR image grldpoint. A mean rainfall rate, 
as specified in the WSR Handbook (NOAA, 1979), was assigned to the entire area 
covered by each radar contour level. The fractional area covered by that 
contour level within each 22.2 km grid box was determined using a computer- 
linked digitizing tablet, which allowed the computation of an area-averaged 
rainfall rate for each grid box (and each corresponding SMMR gridpoint). WSR- 
57 radar from the Daytona Beach station covering a portion of Hurricane David 
is included in Fig. 4 for comparison with the SMMR imagery. 

Referring to Fig. 4, it may be generally noted that the quality of the 
restored images tends to decline with decreasing channel frequency. This 
trend is a consequence of the dramatic decrease in resolution of the original 
SMMR imagery with decreasing frequency. At 10.7 and 6.6 GHz the Case A 
restoration method tends to overenhance large, bright features while adjacent, 
less ‘intense regions are underenhanced. This effect is especially pronounced 
in the Case A 6.6 GHz restoration in the vicinity of the Florida coastline, 
where the bright land mass is abutted by severely underenhanced rainbands. 
Degradation with decreasing channel frequency is also observed in the Case B 
•imagery. In Case B the structure of the rain patterns (at 37 GHz resolution) 
is maintained through the influence of the spatial constraint; however, the 
brightness temperature values at 6.6 GHz show only modest consistency with 
either the radar traces or the imagery in the higher-frequency SMMR channels. 

The poor quality of the 6.6 GHz restorations can be generally atrreibuted 
to the large number of degrees of freedom in the governing system of equation 
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{Eq. 2). This problem was previously discussed in Section 4a. At 6.6 GHz 
there are a large number of degrees of freedom and therefore many possible 
solutions. The restoration procedure iterates to the solution which Is 

closest to the Initial guess. Since the Initial guess field (the original 
blurred image) is quite smooth in Case A, only the large-scale image features 
are enhanced in the final solution. Incorporation of the model constraint 
brightness temperatures in the initial guess of the unconstrained image 
deconvolution did not improve the restored image quality. 

The observed large-scale pattern or over»underenhancement was reproduced 
in the restoration of synthetic images when the size of the blurring function 
template was not small with respect to the image size. The SMMR swath 
dimensions limited the image size to a 6° x 6° grid domain in the present 

study. The simulations also demonstrated that the imposition of accurate 
constraints helped to improve restored image quality. Unfortunately, the 
model constraint of Case B has limited accuracy. Instrument noise or 

calibration errors in the 37 GHz data can propagate through the constraint 

model, producing error in the constraint brightness temperatures at the lower 
frequencies. 

Despite some degradation with decreasing frequency, the restorations at 
37, 21, 18, and even 10.7 GHz show significant improvement in resolution over 
the original imagery. Moreover, there is reasonable consistency between 
channels in the position and intensity of enhanced-resolution features. The 
positions of radar-defined rainbands are well correlated with bright bands in 
these rain-sensitive microwave channels. 

The restoration of Hurricane David also "’fll inter-ch>nnel differences 
which might aid in the retrieval of maritime rai.ifall rates. The theory of 
microwave radiative transfer through extensive precipitating clouds over the 
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ocean Indicates little change In the observed 37 GHz brightness temperature in 
response to changes In rainfall rate at rainfall rates greater than roughly 10 
mm/hr (Olson, et al., 1982). Models also indicate increasing brightness 
temperature sensitivity to rainfall rate as the channel frequency decreases. 
Indeed, evidence of increasing rain sensitivity with dereasing channel 
frequency is seen in the restored imagery. Hurricane eye and rainband 
structure which is not obvious at 37 GHz is revealed at 21, 18, and 10.7 GHz. 

A second example was selected to demonstrate the performance of the 
restoration procedure in the vicinity of a bright land mass. For this 
example, WSR-57 radar and original /restored SMMR imagery of squall line off 
the coast of the Florida panhandle at 1706Z, 7/11/79 are presented in Fig. 
5. The restored SMMR images on the right-hand side of the figure were 
produced using the Case A (unconstrained) and Case B (constrained) enhancement 
methods. 

Although the squall line is evident in the Apalachicola radar imagery, 
blurring of the adjacent Florida coastal region tends to obscure this feature 
in the original 21, 18, 10.7 and 6.6 GHz SMMR data. The problem of signal 
contamination in the vicinity of land masses has contributed to omission of 
such near-coast SMMR data in other studies (Milman, 1983; Wilheit, et al., 
1984). As seen in Fig. 5, the Case A restoration procedure both enhances the 
squall line and distinguishes it from coastal Florida at 37, 21, and 18 GHz. 
Blurring of the bright land mass at 10.7 and 6.6 GHz is so great that 
application of the restoration procedure does not help to delineate coastal 
rain features in the Case A imagery. Application of the model constraint in 
Case B helps to outline rain features at 10.7 and 6.6 GHz, but the brightness 
temperature values in the rain regions, as noted in the preceding example, are 
somewhat suspect. 
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6. CONCLUDING REMARKS 

The nonlinear programming method of Yeh, et al. (1984) performed 
satisfactorily when applied to the restoration of 37, 21 and 18 GHz SMMR 
Images. Utilization of a physically-based constraint helped to delineate rain 
features In restored 10.7 and 6.6 GHz Imagery, but in general the restoration 
of Imagery at these frequencies was difficult, owing primarily to the 
extremely low resolution of the original imagery and the large number of 
degrees of freedom In the image deconvolution problem. 

Constraint effectiveness could be increased if visible, infra-red, or 
other high resolution satellite data were incorporated into the constraint 
area classification shceme. The additional information provided by these 
sensor data could be used to reduce the number of misclassified grid elements 
in the constraint. The accuracy of a constraint based upon the inversion of 
satellite data, such as the one developed in this study, would also benefit 
from future efforts to reduce instrument noise/calibration error. Another 
related subject under current investigation by C.-L. Yeh is the derivation of 
an optimal constraint weighting factor A, which would vary as a function of 
the relatived levels of noise in the original image and the constraint. 

The quality of the restored imagery might also be improved if the antenna 
blurring functions were more accurately estimated. In the present study only 
a limited amount of imagery was used to make blurring function estimates. 
Several coastal swaths in regions where the land protion contains nearly 
homogeneous, dry terrain would be required for accurate estimates. In 
addition a swath-oreinted gridding scheme, rather than the latitide/longitude 
grid employed in this study, would help ameliorate the restoration 
procedure. Blurring functions determined in such a grid system would be 
independent of satellite orbital position. 
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The restoration method presented here may also be applied to other types 
of metorological satellite imagery. Prerequisites for application are (a) the 
data must be interpolated, if necessary, to a regular grid, and (b) an 
estimate of the antenna blurring function must be made In that particular grid 
coordinate system. As it was demonstrated in Section 3, the antenna blurring 
function can be determined using the same nonlinear programming method that is 
utilized to enhance image resolution. In this case one must obtain a subset 
of images for which the ideal -resolution versions can be estimated. This 
approach is especially suited to imagery from multi-resolution sensors such as 
the SMMR or the forthcoming Mission Sensor Microwave Imager (SSM/I), since all 
channels are co-navigated, and there are strong correlations between features 
in different channels. 
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Table 1. Empirical coefficients K v and exponents y v used in Eqs. 10 and 12, 
respectively. The constants y v were derived from Rayleigh theory, assuming 
a homogeneous, Isothermal cloud layer at 273.15 K. 


'[GHZ] 

K v 



37 

.942 

1.00 


21 

.951 

3.47 x 

10“1 

18 

.950 

2.57 x 

io-i 

10.7 

.931 

9.25 x 

10-2 

6.6 

.918 

3.57 x 

10-2 
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Fig. 1. Schematic for grid domains in the image deconvolution probleu. 



(a) 



Fig. 2. Nimbus -7 SMMR 
ascending pass antenna 
blurring functions for the 
(a) 37 GHz 1, (b) 18 GHz 1, 
and (c) 6.6 GHz 1 channels. 
Horizontal displacements are 
measured from the convolution 
template center in km. 
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Fig. 3. Histograms of SMMR image grid elements in raining (solid line) 
and rain-free (dashed line) regions over the Gulf of Mexico at 0540Z, 

7/17/79. Rain data were derived from the WSR-57 radar at Apalachicola, 
Florida. According to the classification method utilized in this study, 
grid elements for which the 37 GHz polarization difference is greater than 
the single standard deviation threshold (54.1 K) are assumed to be rain-free. 
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fig. 4. (continued) 
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Fig. 5. VSR-57 radar, original and restored SMXR imagery of a squall line near 
Apalachicola, Florida at 1706Z, 7/11/79. Case A is unconstrained; Case B is constrained. 
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